Nonlinear resonant absorption of fast magnetoacoustic waves in strongly anisotropic 

and dispersive plasmas 
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The nonlinear theory of driven magnetohydrodynamics (MHD) waves in strongly anisotropic and 
dispersive plasmas, developed for slow resonance by Clack & Ballai [Phys. Plasmas 15(8), 2310 
(2008)] and Alfven resonance by Clack et al. [A&A 494, 317 (2009)], is used to study the weakly 
nonlinear interaction of fast magnetoacoustic (FMA) waves in a one- dimensional planar plasma. 
The magnetic configuration consists of an inhomogeneous magnetic slab sandwiched between two 
regions of semi-infinite homogeneous magnetic plasmas. Laterally driven FMA waves penetrate the 
inhomogeneous slab interacting with the localized slow or Alfven dissipative layer and are partly 
reflected, dissipated and transmitted by this region. The nonlinearity parameter defined by Clack & 
Ballai (2008) is assumed to be small and a regular perturbation method is used to obtain analytical 
solutions in the slow dissipative layer. The effect of dispersion in the slow dissipative layer is to 
further decrease the coefficient of energy absorption, compared to its standard weakly nonlinear 
counterpart, and the generation of higher harmonics in the outgoing wave in addition to the fun- 
damental one. The absorption of external drivers at the Alfven resonance is described within the 
linear MHD with great accuracy. 

PACS numbers: 52.25.Fi; 52.30.Cv; 52.35.-g; 52.35.Bj; 52.35.Mw 



I. INTRODUCTION 

The problem of interacting fast magnetoacoustic 
(FMA) waves with different magnetic structures is not 
only important in the context of astrophysics and solar 
physics, but also in laboratory plasma devices. Space 
and laboratory plasmas are highly non-uniform and dy- 
namical systems and as a consequence they are a natural 
medium for magnetohydrodynamic (MHD) waves. When 
the magnetic plasma configuration is inhomogeneous in 
the transversal direction relative to the ambient magnetic 
field a phenomenon, known as resonant absorption, oc- 
curs (see, e.g., Appert et al. [l[ an d lonson Q). Some 
of the wave energy can be converted into heat in a thin 
layer which embraces the ideal resonant magnetic surface 
when dissipative processes are taken into account. 

In the context of solar physics, the resonant coupling 
of waves was first suggested by Ionson(author?) [3| as 
a possible mechanism for heating coronal loops. Shortly 
after, several studies on the efficiency of resonant ab- 
sorption in the complicated process of coronal heating 
were published by, eja;., Ionson(author?) [2j, Kupe- 
rus et aZ.(author?) [J, Davila(author?) [| and Holl- 
weg(author?) Qj. The same principle was used to ex- 
plain the observed loss of power of acoustic oscillations 
in the vicinity of sunspots by, e.g., Hollweg(author?) 
0, Lou(author?) [8], Sakurai et al. (author?) 
Goossens and Poedts(author?) [T(ij |. Goossens and Holl- 
weg(author?) [11] and Stenuit et aL(author?) [l2l ]. 
All these studies dealt with the Alfven resonant posi- 
tion. Although happening at lower frequencies, slow res- 
onance is also important as shown in a study by Kep- 
pens (author?) [13| where he investigated the interac- 
tion of sound waves with hot evacuated magnetic fibrils. 



Most of the analytical studies of resonant absorption were 
based on the linear theory due to its relative simplicity. 

A new approach to the problem of resonant absorp- 
tion in the context of high Reynolds number plasmas was 
given by Ruderman et al. (author?) [l4| who developed 
a nonlinear theory of resonant absorption for slow waves 
in isotropic plasmas. They pointed out that nonlinearity 
has to be taken into account under typical solar condi- 
tions near resonance. The theory of nonlinear resonant 
slow waves was extended to strongly anisotropic plasmas 
in Ballai et al. (author?) [lj| to describe conditions typ- 
ical for the solar chromosphere and corona. Over the 
next few years there was an enormous amount of effort 
put into studying resonant absorption, including the in- 
vestigation of the effect of equilibria flows at the slow res- 
onance (see, e.g., Ballai and Erdelyi [Hj]), the absorption 
of sound waves at the slow dissipative layers in isotropic 
and anisotropic plasmas (see, e.g., Ruderman et al. [17[ 
and Ballai et al. [HI]) and the effect of an equilibrium 
flow on the absorption of sound and FMA waves due to 
the coupling in the slow continua (see, e.g., Erdelyi and 
Ballai [19j and Erdelyi et a l. l20l| ). In a recent paper, 
Clack and Ballai(author?) [21] showed that in strongly 
anisotropic and dispersive plasmas the dispersion, dissi- 
pation and nonlinearity are all of the same order inside 
the dissipative layer. 

A study by Clack et al. (author?) [22| on the nonlinear 
effects at the Alfven dissipative layer found that nonlin- 
earity and dispersion are always negligible in comparison 
to the linear terms describing dissipation. This implies 
that the linear theory is always applicable for resonant 
absorption at the Alfven resonance if the dimensionlcss 
amplitude of perturbations inside the dissipative layer are 
less than unity. Moreover, Clack et al. (author?) [22j 



showed that the largest nonlinear and dispersive terms 
cancel out - leaving only small corrections to linear the- 
ory. 

Many studies of resonant absorption considered only 
the sound (or slow) and Alfven waves as excellent candi- 
dates for coronal heating. Alfven waves can only carry 
energy along the magnetic field lines and slow waves are 
only able to carry 1 — 2% of energy under coronal (low 
plasma-/3) conditions. However, FMA waves might have 
an important role in explaining the coronal temperatures, 
as has been shown by, e.g., Cadez et al. (author?) (23| 
and Csfk et al. (author?) [241 ]. 

The aim of the present paper is to study the nonlinear 
(linear) resonant interaction of externally driven FMA 
waves with the slow (Alfven) dissipative layer in strongly 
anisotropic and dispersive static plasmas. The governing 
equations and jump conditions derived earlier by Clack 
and Ballai(author?) [U and Clack et aZ.(author?) [H 
will be used to study the efficiency of absorption at the 
slow and Alfven resonance. The paper is organized as 
follows. In the next section we introduce the governing 
equations, the equilibrium state and the fundamental as- 
sumptions which allow analytical progress. In Sec. Ill 
we find the solutions describing the waves outside the 
dissipative layers. Section IV is devoted to the nonlinear 
solution inside the slow dissipative layer. In Sec. V we 
derive the solution inside the Alfven dissipative layer. In 
Sec. VI we will calculate the absorption coefficient in the 
case of slow/Alfven resonance. Finally, in Sec. VII we 
summarize our results and draw our conclusions. 



II. GOVERNING EQUATIONS AND 
ASSUMPTIONS 

The dynamics and absorption of the waves will be stud- 
ied in a Cartesian coordinate system. The equilibrium 
state is shown in Figure [T] The configuration consists 
of an inhomogeneous magnetized plasma < x < xq 
(Region II) sandwiched between two semi-infinite homo- 
geneous magnetized plasmas x < and x > xq (Regions 
I and III, respectively). We have chosen this model to ob- 
tain analytical results. Our intention is to have a model 
which gives us the trend in the absorption of an inci- 
dent wave on a magnetic structure. It is obvious that 
real magnetic structures are more complicated (and far 
from being fully understood), however, the magnetic field 
has been simplified to be unidirectional in order to make 
the model more transparent, such that the role of the 
dispersion at the resonance and the change in the ab- 
sorption can be investigated more fully, and compared 
to previous studies. We took inspiration for this model 
from seminal studies such as Ruderman et al. (author?) 
[13], Ballai et al. (author?) [18], Erdelyi(author?) [H, 
Roberts(author?) [25|, Edwin and Roberts(author?) 
[H] and Ruderman(author?) (27[. 

The equilibrium density and pressure are denoted by 
p and p. The equilibrium magnetic field, B, is unidi- 
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FIG. 1: Illustration of the equilibrium state. Regions I (x < 0) 
and III (x > 0) contain a homogeneous magnetized plasma 
and Region II (0 < x < xo) an inhomogeneous magnetized 
plasma. The shaded strip shows the dissipative layer embrac- 
ing the ideal resonant position x c . 



rectional and lies in the yz-planc. In what follows the 
subscripts "e" , "0" and "i" denote the equilibrium quan- 
tities in the three regions (Regions I, II, III, respectively). 
It is convenient to introduce the angle, a, between the z- 
axis and the direction of the equilibrium magnetic field, 
so that the components of the equilibrium magnetic field 
are; B y = _Bsinaandi? 2 = Bcosa. All equilibrium 
quantities are continuous at the boundaries of Region II, 
so they satisfy the equation of total pressure balance. It 
follows from the equation of total pressure that the den- 
sity ratio between Regions I and III satisfy the relation 



Pi 

Pe 



24, 



(1) 



where the squares of the Alfven and sound speed are v A = 
Bq/hqPo and c| = jpa/pa- Where po is the magnetic 
permeability of free space and 7 is the adiabatic constant. 
Replace the subscript "0" with "e" for Region I and a i" 
for Region III. We consider a hot magnetized plasma such 
that c 2 Si > cf e , and v% > v\ e . 

The objective of the present paper is to study (i) the 
combined effect of nonlinearity and dispersion on the in- 
teraction of incoming fast waves with slow dissipative lay- 
ers and (ii) the interaction of incoming fast waves with 
Alfven dissipative layers. We, therefore, have two dif- 
ferent criteria. For interaction of FMA waves with the 
slow dissipative layer we assume that the frequency of the 
incoming fast wave is within the slow continuum of the 
inhomogeneous plasma, so that there is a slow resonant 
position at x — x c in Region II. Interactions with the 
Alfven dissipative layer leads to the assumption that the 
frequency of the incoming fast wave is within the Alfven 
continuum of the inhomogeneous plasma, so that there 
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is an Alfven resonant point at x = x a in Region II. This 
leads to the inequality, c-r e < Lo/k < en , at the slow reso- 
nance. Where the square of the cusp speed, c\, is defined 



dispersion relation ([3]) becomes 



by 



4«i/(4 + v 2 A ) 



We also have the inequality, 



VAe < uj/k < i>Ai, for Alfven resonance. Here u> is the fre- 
quency of the incoming fast wave and k = (k 2 + k 2 ) 1 ^ 2 is 
the wave number. Even though, in principle, when a slow 
resonance occurs in this manner an Alfven resonance is 
also present we ignore the Alfven resonance that occurs 
alongside the slow resonance as this would complicate the 
analysis and obscure the results associated with the slow 
resonance. We study the Alfven resonance separately to 
the slow resonance. We note that the Alfven resonance 
would, in simple terms, act to restrict the energy avail- 
able at the slow resonance. We intend to address the issue 
of coupled resonances in our next paper, where we will 
show that the governing equations derived here remain 
the same (meaning the work here is valid), however, the 
interaction of the waves between the resonant positions 
changes the absorption of wave energy. 

In an attempt to remove other effects from the analysis 
we consider the incoming fast wave to be entirely in the 
xz- plane, i.e. k y — 0. Ruderman et al. (author?) [l7| 
suggests aligning the equilibrium magnetic field with the 
z-axis, to remove the Alfven resonance (if we consider 
planar waves) from the analysis for slow resonance, how- 
ever, this is not possible nor necessary here. The disper- 
sion is dependent on the angle between the equilibrium 
magnetic field and the z-axis (a), hence if a = the dis- 
persion effects disappear, and we recover t he g overning 
equation studied by Ballai et al. (author?) (18j ]. 

The inequalities above guarantee that the slow and 
Alfven resonances appears in Region II when studying 
in the upper chromosphere and the solar corona, re- 
spectively. The resonant positions, therefore, are de- 
fined mathematically as: lo c = kcT(x c ) cos a and uj a = 
kvA(x a ) cos a. The position of the resonant points also 
provides us with some information about the plasma con- 
dition. First, in conjunction with Eq. |T]) we obtain that 
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Hence, the plasma in region III is more rarefied than in 
Region I. Secondly, it follows that cx e < c n and the 
plasma in Region III is hotter than the plasma in Region 
I. 

The dispersion relation for the impinging propagating 
fast waves takes the form 



k 2 



4) 



4) 2 



iv A c 2 s cos2 1 



1/2 



(3) 

where <fi is the angle between the direction of propagation 
and the background magnetic field within the rcz-plane 
and k = k x e x + k z e z . For the sake of simplicity, we 
denote the ratio k x /k z . Since the equilibrium mag- 
netic field in the zz-plane is aligned with the z-axis, the 
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where 1 + n 2 = 1/ cos 2 <f>. 

We assume the plasma is strongly magnetized in the 
three regions, such that the conditions Wj( e )Tj( e ) 3> 1 are 
satisfied, here Wj( e ) is the ion (electron) gyrofrequency 
and 7"j( e ) is the ion (electron) collision time. Due to 
the strong magnetic field, transport processes are derived 
from Braginskii's stress tensor (see, e.g., Braginskii|28|: 
Ruderman et al.^K^). As we deal with two separate 
waves (slow and Alfven), we will need to choose the 
particular dissipative process which is most efficient for 
these waves. For slow waves, it is a good approximation 
to retain only the first term of Braginksii's expression 
for viscosity, namely compressional viscosity(author?) 
[30| . In addition, in the solar upper atmosphere slow 
waves are sensitive to thermal conduction. In a strongly 
magnetized plasma, the thermal conductivity parallel to 
the magnetic field lines dwarfs the perpendicular com- 
ponent, hence the heat flux can be approximated by the 
parallel component only(author?) [3l[. On the other 
hand, since Alfven waves are transversal and incompress- 
ible they are affected by the second and third compo- 
nents of Braginskii's stress tensor, called shear viscos- 
ity (author?) [22j. Finally, Alfven waves are efficiently 
damped by finite electrical conductivity, which becomes 
anisotropic under coronal conditions. The parallel and 
perpendicular components, however, only differ by a fac- 
tor of 2, so we will only consider one of them without 
loss of generality. All other transport mechanisms can 
be neglected. For further details, please refer to, for ex- 
ample, Clack et aL(author?) |22]|. BraginskiKauthor?') 

28 1, Ruderman et al. (author?) [29], Hollweg(author?) 

3C|, Priest (author?) [3l[ and Porter et al. (author?) 

The dynamics of nonlinear resonant MHD waves in 
anisotropic and dispersive plasmas was studied by Clack 
and Ballai(author?) [2l| and Clack et al. (author?) 
[22| . They derived the governing equations and con- 
nection formulae necessary to study resonant absorp- 
tion in slow/Alfven dissipative layers. We recall the 
key steps and necessary results found by Clack and Bal- 
lai(author?) [H| and Clack et al. (author?) [12]. 

The efficiency of dissipation, when studying slow dissi- 
pative layers, in an anisotropic plasma is given by the 
(compressional) viscous Reynolds number {R e / C )) and 
the Pechlet number (P e ), combining to define the to- 
tal Reynolds number: R^ 1 = R~^ c \ + P^ 1 , where i? e ( c ) 
and P e are defined by i? e ( c ) = v c hl cn po c /Vo and P e = 

v c hlchPocR/H\\- Here v cn is the characteristic velocity 
(e.g. the slow magnetoacoustic velocity at x = x c ), l cn 
is the characteristic length, p 0c = p (x c ), R denotes the 
gas constant and ku is the coefficient of thermal con- 
ductivity parallel to the equilibrium magnetic field lines. 
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The efficiency of dissipation, when studying Alfven dis- 
sipative layers, in an anisotropic plasma is measured in 
a slightly different way. Now dissipative processes are 
described by the (shear) viscous Reynolds number (R e ) 
and the magnetic Reynolds number (R m ), combining to 
define the total Reynolds number: R~ x = R^ 1 + i?" 1 , 
where R e = v ch l ch po a /m and Rm = v ch l ch /fj. Here v ch 
is the characteristic velocity (e.g. the Alfven velocity 
at x = x a ), l c h is the characteristic length, r\\ is the 
coefficient of shear viscosity and ij is the coefficient of fi- 
nite electrical resistivity. Originally, these total Reynolds 
numbers were introduced based on intuition, simplicity 

and linear theory (see, e.g., Sakurai et al. I9(, Goossens 

et al. |34[ and Goossens and Ruderman [35|]). However, 
it turned out that using these definitions the strength of 
dissipation is the same order of magnitude as the inverse 
of the total Reynolds numbers. Under chromospheric and 
coronal conditions R 3> 1 which means that dissipation 
is only important inside the dissipative layer. Far away 
from the dissipative layer amplitudes are small, therefore 
we can use the linear ideal MHD equations to describe 
the plasma motions far from the resonant position. These 
equations can be reduced to a system of coupled first or- 
der PDE's for the total pressure perturbation, P, and the 
normal component of the velocity, u, 



du V dP BP 
dx F d6 ' dx 



PqA du 



(5) 



Here 
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:) (V 2 - 
A = V 2 - v\ cos 2 a. 
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The system ([5]) describes the wave motion far from the 
ideal resonant position. The singularities in the coeffi- 
cients A and C give the conditions of Alfven and slow 
resonance. All perturbations depend on the combination 
9 = z — Vt, where V = oj/k is the phase speed. 

Inside the thin dissipative layers (where the dynamics 
is described by the nonlinear and dissipative MHD equa- 
tions) embracing the ideal resonant surfaces (x = x c , x = 
x a ) we must use the governing equations derived by Clack 
and Ballai(author?) [H| and Clack et aZ.(author?) [13]. 
The characteristic thickness of the slow dissipative layer, 
5 C , is 



5r = 



V 3 k\ 



Ac|(c; 



Sc 



(7) 
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Here k = 2ir/L with L the wavelength, the subscript "c" 
indicates that the quantity has been calculated at the 
slow resonant position. The quantity A is defined by 



7?o(2^ c +3c| c ) ( 7 - 1) 2 K|| [y\ c + 4 e 



jpo c Rcl r 



(8) 



and A c is simply the gradient of the cusp speed given by 
A c = — (dc^/dx) cos 2 a. Clack and Ballai(author?) 
PH showed that nonlinearity and dispersion are impor- 
tant in the slow dissipative layer if the nonlinearity pa- 
rameter is greater than unity, N 2 = eR 2 > 1, where e 
is the dimensionless wave amplitude far from the dissi- 
pative layer. The concept of nonlinear parameters was 
introduced by Ruderman et al. (author?) [14| for slow 
waves and Clack et aZ.(author?) [22| for Alfven waves. 
The two parameters are different not only in their form 
but also in the values the Reynolds numbers take. In the 
case of slow waves (damped by compressional viscosity, 
i.e. the first term in the Braginskii's viscosity tensor) 
the Reynolds number that corresponds to a characteris- 
tic length of 200Mm, a speed of 200kms _1 , a density of 
10 -13 kgm -3 and a compressional viscosity coefficient of 
5 x 10 _2 kgm _1 s _1 is about 80. Alfven waves are effi- 
ciently damped by shear viscosity which is given by the 
second and third coefficients of the Bragisnkii's tensor 
(here denoted cumulatively as rji). Since r\\ = r]o / \ujiTi) 2 
and under coronal conditions uiiTi is of the order of 10 5 , 
we obtain that the coefficient of shear viscosity is about 
10 orders of magnitude smaller than the coefficient of 
compressional viscosity. Now, using the characteristic 
speed of lOOOkms -1 , the Reynolds number used in cal- 
culating the nonlinear parameter in the case of Alfven 
nonlinearity is 4 x 10 12 . The nonlinearity parameter for 
resonant Alfven waves is eR 2 ' 3 <C 1. However, it was 
shown by Clack et al. (author?) [22| that the waves in 
this situation remain linear anyway. In the present pa- 
per, therefore, we do not need the nonlinearity parameter 
for resonant Alfven waves. The characteristic thickness 
of the Alfven dissipative layer, 5 a , is 



5 a 



V 



k\A a \ 



POa 



1/3 



(9) 



with A a being the gradient of the Alfven speed given by 
A a = — (dv 2 A /dx) a cos 2 a. 

The governing equation inside the slow dissipative 
layer is(author?) [2l[ 



dq c , . dq c ._ 1 d 2 q c 

Gc ^e + Aqc ^e - k -cW 

dq c dq c 



kV 4 dP 



da 89 p 0c v 2 Ac \A c \ dO 



where 



A = R 



5 C ' 

2 v\ JA\ [(7+1K.+34,] 
kV 8 



* = R 2 



x\^\ 4 e v A c ( v A c + c D sin< 



kV 13 



X = -quj e Tc. 



, (io) 

(11) 

(12) 

(13) 
(14) 



Here the first term of the governing equation appears 
due to the inhomogeneity in the cusp speed, the second 
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term describes the nonlinearity of waves, the third term 
stands for the dissipative effects while the last term on 
the left-hand side describes the nonlinear dispersive ef- 
fects generated after taking into account Hall currents by 
Clack and Ballai(author?) [2l|. The term on the right- 
hand side can be considered as a driver. We also note 
that q c (<7 c , 9) is the dimensionless component of velocity 
parallel to the equilibrium magnetic held and y = rjuj e T e 
is the coefficient of Hall conduction(author?) [2lj . 

The governing equation inside the Alfven dissipative 
layer is(author?) [22| 

dq a _ , d 2 q a _ A; sin a dP 
aa d8 k dal -p 0a \A a \de> 

with a a = {x—x a )/8 a . Here q a (a a ,9) is the dimensionless 
component of velocity perpendicular to the equilibrium 
magnetic held. We should point out here that although 
nonlinearity and dispersion have been considered when 
deriving the dynamics of the Alfven resonance, the gov- 
erning equation remains linear regardless of the degree of 
nonlinearity (for details see Clack et al. (22J). 

When studying resonant MHD waves, we are generally 
not interested in the solution inside the dissipative layer 
and can consider the dissipative layer as a surface of dis- 
continuity. Instead, we solve the system ([5]) and match 
the solutions at the boundaries of the discontinuity using 
connection formulae. These connection formulae deter- 
mine the jumps in u and P across the dissipative layer. In 
the context of solar plasmas, they were first introduced by 
Sakurai et as?. (author?) [§]. It was shown by Clack and 
Ballai(author?) [2l[ (in complete agreement with Rud- 
erman et a?, (author?) [TJ] and Ballai et a?, (author?) 
[ll|) that the first connection formula is [P] — 0, where 
the square brackets denote the jump across the dissipa- 
tive layer. It can also be shown, in a similar manner, that 
the same jump condition exists for the Alfven resonance. 
The second connection formula for slow resonance can 
only be written in implicit form, i.e. 

m = ~~u~~~2 — ^ r % da > 

fccos^ a J_ oc 09 

where we use the Cauchy principal value of the integral 
because the integral is divergent at infinity. As a result 
we must solve Eqs. <j5j> and (jTUJ) along with the boundary 
conditions, [P] = and Eq. (jTHJ) . In an attempt to fol- 
low the same procedure utilized for finding solutions at 
the slow resonance we can write the jump in the normal 
component of velocity for the Alfven resonance in an im- 
plicit form. For the sake of brevity, we do not show the 
derivation here, but it follows the procedure to find the 
jump in the normal component of velocity completed by 
Clack and Ballai(author?) [21]. This jump is given by 

Finally, we should note some critical assumption we 
make to allow analytical progress. From the very begin- 
ning we must assume that the nonlinearity parameter is 



small so that regular perturbation theory can be applied 
at the slow resonance. We also assume that the inho- 
mogeneous region is thin in comparison with the wave- 
length of the impinging wave, i.e. kxo <C 1. Ruder- 
man(author?) [27| investigated the absorption of sound 
waves at the slow dissipative layer in the limit of strong 
nonlinearity. In his analysis nonlinearity dominated dis- 
sipation in the resonant layer which embraces the dis- 
sipative layer. He concluded that nonlinearity decreases 
absorption in the long wavelength approximation, but in- 
creases it at intermediate values of kxo, however, the in- 
crease is never more than 20%. To the best of our knowl- 
edge, at present, we cannot solve the governing equation 
(|10p in the limit of strong nonlinearity due to the nonlin- 
ear dispersive term, therefore we restrict our analysis to 
the weak nonlinear limit. We mention that no such as- 
sumptions are needed for studying the Alfven dissipative 
layer since the governing equation (|15[) is linear. 



III. SOLUTIONS OUTSIDE THE DISSIPATIVE 
LAYERS 

In what follows we derive a solution for the system 
([5]) in Regions I, II and III. In Region II we only find 
the solution outside the dissipative layers. Section IV is 
devoted to finding a solution to Eq. (JTUJ) inside the slow 
dissipative layer and Section V is used to hnd a solution 
to Eq. (Tl5|) inside the Alfven dissipative layer. Outside 
the dissipative layers, the solutions take identical forms. 



A. Region I 

The solution of Eq. ((5|) in Region I is given in the form 
of an incoming and outgoing fast wave of the form 

P = e{p e cos [k (9 + K e x)\ + Acqs [k{9 - K e x)]} , (18) 
K e V {p e cos [k (9 + K e x)] — A cos [A; (9 — «s e a;)]} 
Pe (V 2 - v 2 Ae cos 2 a) 

(19) 

where e <C 1 is the dimensionless amplitude of perturba- 
tion far from the dissipative layer. The frequency of the 
incoming wave is given by Eq. (|H) and must lie within 
the slow or Alfven continuum depending on which dissi- 
pative layer we are studying. The first term in Eqs. (fTS"]) 
and (|19p describes the incoming wave, while the second 
term describes the outgoing wave which will be obtained 
in Section IV for slow dissipative layers and in Section V 
for Alfven dissipative layers. 
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B. Region II 

In Region II, the equation for the total pressure, P, is 
obtained by eliminating u from the system ([5]), 



dx 



dP 



Po (V 2 — v\ cos 2 a) dx 



8 2 P 



(20) 



Since we have assumed kxa <C 1, the ratio of the right- 
hand side and the left-hand side is of the order of k Xq. 
It follows that 

dP 

— = po {V 2 - v\ cos 2 a) f(9) + 0{k 2 x 2 ), (21) 

where the function f(9) is determined by the second 
equation of ^ and the boundary conditions at s = 0. 
Equation (|2"Tj) yields 

P = P{6) + f(9) [ po [V 2 - v 2 A cos 2 a] dx + 0{k 2 xl). 
J o 

(22) 

The function P(9) has to be determined by the bound- 
ary conditions at x = 0. It can be shown that, because 
[P] = 0, the functions f(9) and P{9) take the same val- 
ues throughout Region II. Noting that the second term 
in Eq. ([22"]) is of the order of kxo we can express P in a 
simplified form 



P = P{9) + (kxo) P'{x, 9) + 0{k 2 xl). 



(23) 



IV. WEAK NONLINEAR SOLUTION INSIDE 
THE SLOW DISSIPATIVE LAYER 

Since we are not able to solve the governing equation 
(fTT)|) inside the slow dissipative layer analytically, we con- 
sider the limit of weak nonlinearity (N 2 <C 1). In ac- 
cordance with this assumption we rewrite the governing 
equation (JTDJ) and the jump condition p6[) as 



dqj -i7A\_(^ -x^9q^dq c 
a ^9 +£ C Uh c W £ C ^^9 



-id 2 q, 



V 4 



dP r 



89 2 



p 0c wl c |A c |a;o dO 



vc °^r?±d*, 



, (26) 



(27) 



where 



Tc = ^r o , C=^g^, D 2 = eR 4 = R 2 N 2 . (28) 

Note that £ is of the order of eR 2 , the ratio (A /if?) is of 
the order of unity and qZ is of the order of e. In what 
follows we drop the bar notation and for the rest of this 
section we drop the subscript "c" on the dimensionlcss 
variable q. 

We proceed by using a regular perturbation method 
and look for solutions in the form 



C. Region III 

To derive the governing equation for Region III we 
eliminate the normal component of the velocity from the 
system ([5]) to arrive at 



8 2 P 
dx 2 

where k 2 is defined as 



+ K? 



, d 2 P 

d9 2 



= 0, 



V* V 2 (c 2 Sl + v 2 Al ) + , 



'■Si V Ai 



cos a 



(4i + <)(v 2 



■Ti 



cos 2 a) 



(24) 



(25) 



Since, for slow dissipative layers, V < CTi cos a, it follows 
that k 2 > 0. It also follows that for Alfven dissipative 
layers k 2 > because V > ct% cos a > vac cos a. There- 
fore, Eq. (HM)) is an elliptical differential equation and 
the wave motion is evanescent in Region III. In reality, 
there could be wave leakage. The existence of wave leak- 
age depends on the profile of the slow and Alfven speeds 
in the inhomogeneous region (Region II). For simplicity, 
we have assumed that the slow and Alfven resonances 
take place at a single location (obviously different for the 
two resonances), which means the profiles of the slow and 
Alfven speeds are monotonically increasing inside Region 
II. Should we have a more complex model, the possibility 
of wave leakage would need to be taken into account. 



f = eY / C- 1 fn, (29) 



where / represents any of the quantities P, u and q. 



A. First order approximation 

In the first order approximation, from Eq. (|26|). we 
obtain 



9qi , _i d 2 q 1 
cr-r— k 



V 4 dP lc 



89 " 89 2 ~ P ocV Ac \k\xo d9 



(30) 



Since the total pressure, P, is continuous throughout the 
dissipative layer and is periodical with respect to 9, we 
look for a solution in the form gi — $l(gie lke ), where g\ 
represents Pi, Ui and q\ and 5R indicates the real part of 
a quantity. 

In Region I the solutions for the pressure and velocity 
exactly recover the results found in linear theory, i.e. 



i = p e e + A\e 



K e V (p e e ikK - x - A ie - lkK ' x ) 
p e (V 2 - v 2 Ae cos 2 a) 



(31) 
(32) 



where Ai (and subsequent values of Aj/j is the amplitude 
of the outgoing wave. The first terms of the right-hand 
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side of Pi and Ui represent the incoming wave, while 
the second terms are the outgoing (reflected) wave. The 
continuity of the total pressure perturbation at x = 
and x = xo in combination with Eq. (|23p yields Pi, in 
Region II, as 



Pi =p e +A x + (fcx )/ii, 



(33) 



where h n = h n {x) = P^(x) - P^(Q), n > 1. The so- 
lution in Region III is obtained by using Eqs. (JSJ, (|24]1 
and (|33p with the continuity conditions at a; = xq. The 
solution takes the form 



"1 



Pi (U 2 - tr^ cos 2 a) 



(34) 



Utilizing the fact that u\ is continuous at x = and 
x = xo, and employing Eqs. §5§ and (|33[) we find that 
the jump in the normal component of velocity across the 
dissipative layer is 



iKiV(p e + Ai) 



K e V(p e - Ai) 



Pi (V 2 - v 2 Ai cos 2 a) p e (V 2 - v 2 Ae cos 2 a) 
- ikV (p e + Ax) W F- 1 {x)dx 



«i 



r x ° hi (x) 

-ikV(kx )& / dx ' ( 36 ) 

Jo F(x) 

where the expression of F(x) is given by Eq. ([5]). 
Solving Eq. ([30]) reveals qi to be 



U 4 (p e + Ai){l + ^(fca;o)} 



(37) 



Substitution of this result into Eq. I|27p leads to another 
definition of the jump in the normal component of veloc- 
ity across the dissipative layer, namely, 



-TTkV 5 (j> e + A X ) {1 + 0{kx O )} 

p ac v Ac \A\ cos 2 a 
Comparing Eqs. (f3l>|) and (l3"5)) we obtain that 

a t - p + iv ^,,2 2^ 



T + p + IV 



(38) 



(39) 



where 



tt/cU 5 



p 0c ^ c |A|cos 2 a' p e (V 2 - v 2 Ae cos 2 a 



Pi (V 2 - v 2 Ai cos 2 a) 



W3? I F^ix) dx. (40) 



When deriving Eq. (I39p we have employed the estimate 
that k&> f*° h n (x)/F{x) dx = 0(kx o ). The quantity A x 
is a complex value. This means that the outgoing (re- 
flected) wave has a phase alteration compared with the 



incoming wave. The true amplitude of the outgoing wave 
is given by A\ = {A\,,+A 2 , im y) x / 2 (where the subscripts 
"r" and "im" mean the real and imaginary parts, respec- 
tively). The Fourier analysis allows Ai to be complex. In 
general, a complex value of A n means the true amplitude 
of the outgoing harmonic is defined as above and a phase 
of the outgoing wave is shifted by tan -1 ^^/^^). 
This definition of A n applies to all subsequent orders of 
approximation. 

In Ruderman et al. [TtJ and Ballai et al. [H| a similar 
procedure was carried out. Our results are similar with 
theirs if we consider B e = and a = 0. This conclusion 
is not surprising because the first order approximation 
with respect to the nonlinearity parameter coincides with 
linear theory. In addition, dispersion due to the Hall 
effect at the slow resonance does not alter linear theory 
either since dispersion effects appear as a nonlinear term 
in the governing equation. 



B. Second order approximation 

Nonlinear effects start to be important from the second 
order approximation onwards, but they are always due to 
the nonlinear combination of lower order harmonics. In 
this order of approximation Eq. (|26[) is reduced to 



dq 2 ,-i<9 2 g 2 



u 4 dP 2 , 



d0 2 p 0c v Ac \k\x dd 



dqx dqi dqi 

+ (41) 



Taking advantage of the form of the first order approx- 
imation terms enables us to rewrite the second term on 
the right-hand side of this equation as 



dqi „ / ik . 2 2ike \ 



(42) 



Since the nonlinear terms are proportional to K (e 2lfee ) 
it is appropriate to seek a solution of the form g 2 = 
3? (g 2 e 2lke ), where g 2 represents P 2 , u 2 and q 2 . 

Using the same techniques as in the first order ap- 
proximation, it is straightforward to find the jump in the 
normal component of velocity in Region II 



"2 



iKiVA 2 



n e VA 2 



Pi {V 2 - v 2 Al cos 2 a) p e (V 2 - v 2 Ae cos 2 a) 

- 2ikVA 2 0> / F^ix) dx 
Jo 

-2ikV(kx Q )0> -4-T dx. (43) 
Jo F (x) 

Using Eqs. ((37]) and (gSJ we can solve Eq.flHJ to ob- 
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tain 



92 



2 1 



V 4 A 2 



Pocv\ c \A\x 



7 8 (p e + A 1 ) 2 (l + 4fi 2 ) 



Ap 2 c v_ 



Ad 



if 



(44) 



where ft 2 — \/{a — i) is the additional factor due to 
the nonlinear dispersion (as are all subsequent values of 
Oi, i > 2). We substitute the expression for q 2 into Eq. 
PTf to find 



[U2 



2nkV 5 A 2 



P0cV Ac \ 



(45) 



where the terms of the order of k 2 x\ are not indicated. 
To calculate A 2 we compare the jump in the normal com- 
ponent of velocity across the dissipative layer defined by 
Eqs. (O and gg). This leads to A 2 = ff{k 2 x 2 ). This 
result implies that all quantities in the second order ap- 
proximation are zero outside the dissipative layer up to 
an accuracy of 0(kxo). With this restriction the outgoing 
wave remains monochromatic in the second order approx- 
imation. This result coincides with the results of Rud- 
erman et al. (author?) \vh . Ballai et al. (author?) [l8j |. 
Erdelyi et al. (author?) [2(| and Ruderman(author?) 
[27l | (this is especially surprising because in this paper 
nonlinearity is strong). 



C. Third order approximation 

The third order approximation with respect to £ is 
governed by 



dq3 



,_i5 2 Q 3 



V 4 



dP- 



3r 



09 89 2 P0c v 4 Ac \A\x Q d6 

d (gig2) dqi dq 2 dq 2 Oqx 

00 Oa 00 Oa 09 ' 1 ' 

Taking into account the form of the solutions in the pre- 
vious two orders of approximation we can rewrite the 
second term on the right-hand side of Eq. ([4U|) as 



- -M (3iqxq 2 ( 



3ik0 



00 



tq 1 q 2 e 



ik0\ 



(47) 



where q n — 5i (q n e + q*e~ mk0 ) and the asterisk de- 
notes a complex conjugate. This result inspires us to 
seek solutions in the third order approximation in the 
form g 3 = $1 (g3ie lkS + g 33 e 3M ) , where 33 represents 
P3, U3 and (73. Considering the length of this paper we 
only calculate the §31 quantities, as it can be shown that 

A33 = &(k 2 X 2 ). 

In a similar manner as the first and second order ap- 
proximations, we find that the jump in the normal com- 



ponent of velocity across the slow dissipative layer to be 
iKiVA 31 K e VA 31 



[U31 = 



pi (V 2 — v 2 Ai cos 2 a) p e (V 2 — v 2 Ae cos 2 a) 



-ikVA 3 i&> / F-^x) dx 
Jo 





r° (x) 

-ikV(kx )& / -f^-dx, (48) 
Jo F (x) 

To find g 3 i we must exploit Eqs. ([37]), (JUJ) and dUJ 
to solve Eq. (|4"rj)) . The calculation is analogous to the 
first and second order approximation calculations and we 
arrive at the solution 



V 4 A; 



931 



31 



PocV Ac \A\x Q (a 

rl2 



V 12 (p e + A 1 )\p e + A 1 \ 2 (l + 2Q 3 i) 



Sp&ZWtf (a- l y(a-2 l )(a 2 + l) 
where f^i and is given by 

- 3 - '8 + 7i)a 2 - (11 + \2i)a - (44 - hi) 



(49) 



^31 — ^2 



(a - 2i) {a 2 + 1) 



We substitute this expression for 531 to find a second def- 
inition for the jump in the normal component of velocity 
across the slow dissipative layer (up to an accuracy of 
kx ) 



K31J 



nkV 5 A 31 



+ 



P0cV A 

nkV 13 (p e + Ax) \ Pe + A x \ 2 (27 - 8i) 
96pQ C v A 2 c \A\ 3 xl cos 2 a 



, (50) 



Similar to the first two orders of approximation, we can 
compare Eqs. (|4"8|) with ([50)) to find the coefficients A31 



A31 = 



3 - 3 p 3 (27 



Pe T 



8i) cos 4 a 



12w 2 V 2 k 2 x 2 (p + iv) 2 (p 2 + v 2 ) ' 



(51) 



When calculating A31 we have used the estimates r = 
0{kx o ) and kV@> ^° F~ 1 {x) dx = ff{kx ), and retain 
only the terms of lowest order with respect to kxo, as 
we have assumed that kxo <C 1. Equation (I5ip illus- 
trates that with an accuracy of up to €?(kxo) the outgo- 
ing (reflected) wave remains monochromatic in the third 
order approximation. Nevertheless, there is a slight al- 
teration to the amplitude of the fundamental harmonic 
of the outgoing wave from Ax to Ai + C 2 ^3i- These 
results coincide, qualitatively, with the findings by Rud- 
erman et aZ.(author?) [ijj, Ballai et aL(author?) [HI ] 
and Erdelyi et al. (author?) [20], however, A31 is quan- 
titatively larger than that of previous studies and has an 
imaginary component. This implies that the amplitude 
of the wave is greater and the phase of the correction 
is changed when compared with those studies. The ex- 
pression for A 3 i, Eq. (p)Tj) , is different to the ones they 
obtained because of the inclusion of dispersion through 
the Hall current. 
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D. Higher order approximations 



151) , we are left with 



In the fourth order of approximation the outgoing (re- 
flected) wave becomes non-monochromatic. This means 
the energy from this order of approximation no longer 
contribute to the fundamental harmonic, but to a higher 
one. For full details of the calculation please refer to the 
Appendix. 

Continuing calculations to even higher order approxi- 
mations it can be shown that the higher order harmon- 
ics (third, fourth, etc.) are generated in the outgoing 
(reflected) fast wave. The pressure perturbation of the 
outgoing wave can be written as 



P' = eK 




(52) 



The second harmonic only appears in the outgoing wave 
in the fourth order approximation, whereas, higher har- 
monics appear in higher orders of approximation. This 
implies that the estimate A n = <??(C 3 ), n > 2 is valid. 



V. SOLUTION INSIDE THE ALFVEN 
DISSIPATIVE LAYER 

We can find the jump in the normal component of ve- 
locity at the Alfven resonance explicitly, however, in an 
attempt to follow the procedure in the last section (and 
to verify the theory), we proceed to use the implicit form 
of the jump conditions. As the governing equation l|15p 
is linear we only need to calculate one order of approxi- 
mation. 

Although the Alfven resonant position is at x = x a , 
compared with x = x c for the slow resonant position, 
we can use some of the same formulae as in the previ- 
ous section. First, we look for a solution in the form of 
gi = ^.{gie lke ). In Region I, we use Eqs. (01]) and J32]) to 
represent the pressure and normal component of velocity 
perturbations, respectively. For Region II, due to the first 
connection formula, [P] = 0, we can write the pressure 
perturbation as Eq. We also find that Eqs. (131| 

and (|35[) can be used to represent the pressure and nor- 
mal component of velocity perturbations, respectively, in 
Region III. The fact we can employ the same equations 
(as in slow resonance) in the three regions leads to one 
of the definitions of the jump in the normal component 
of velocity over the Alfven dissipative layer being defined 
as Eq. (|36p. It should come as no surprise that this 
definition of the jump across the Alfven dissipative layer 
coincides with the jump across the slow dissipative layer 
in the first order approximation. We are using linear the- 
ory to obtain both expressions and are not looking inside 
the, respective, dissipative layers', so the forms should be 
identical. 

To find q a , so that we find the other definition of the 
jump in Ua-, requires a different approach to the one uti- 
lized in the section before. After Fourier analyzing Eq. 



d 2 q a ifcsina 
da 2 a p a|A Q | 



Pa- 



(53) 



To solve Eq. ([53"]) we introduce the Fourier transform 
with respect to a: 



(54) 



Then from Eq. (153")) we have 



dr A)a|A Q 



. 2-Kiksma (p e + A) ., , , . 
^[la] = nrn ~ S (r), (55) 



where S(r) is the delta-function. We find that the solu- 
tion to Eq. ([55)1 that is bounded for \r\ — > oo is 

2ink sin a(p e + A) 3/3 
&[q a \ = rr—, H{-r)e 1 . (56) 

P0a|A Q | 

Here H(r) denotes the Heavyside function. It was shown 
by Ruderman and Goossens(author?) [17| that 

&> J™ f(a) da= l - ( r Hm o + r fim &[f\\ . (57) 



With the aid of Eqs. U7J), (O and §5T§ we find that 
TrkV(p e + A) sin 2 a 



P0a|A a | 

Comparing Eqs. (|3"6"|) and ([55]) we derive that 
A = - Pe p + 0(k 2 xl), 

T a + + 



(58) 



(59) 



where r a = wkV sin 2 a/ (po a \A a \), and \i and v have their 
forms given by Eq. (|40[) . However, their values are dif- 
ferent for the two resonances. 



VI. COEFFICIENT OF WAVE ABSORPTION 

The coefficient of wave absorption is defined as T — 
(Hin — n out )/Hi n , where n; n and n out are the normal 
components of the energy fluxes, averaged over a period, 
of the incoming and outgoing waves, respectively. It is 
straightforward to obtain that 



r = i-l£KI 2 



C 2 F 



ND, 



(60) 



where Tl is the linear coefficient of wave absorption and 
Tnd is the nonlinear and dispersive correction. Note that 
Tne) is multiplied by the small factor £ 2 which means that 
this term will provide small corrections to linear results. 
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Carrying out calculations we find at the slow reso- 
nance, in agreement with linear theory, that 



r L = 



4t^ 



&{k 2 x%). 



(61) 



The coefficient Tnd is defined as Tnd = 
— (2/j3g)3? {^^31}, which can be rewritten using 
Eqs. (Qll) and ([HI]) as 



ND 



27p 2 e r 6 fx 6 cos 4 



6ir 2 V 2 k 2 XQ (n 2 



+ 0{k 2 x 2 ). (62) 



Both Tl and Tnd are of the order of kxo. This result 
is qualitatively the same as Ruderman et al. (author?) 
[l7l ] and Ballai et aZ.(author?) [l8| results, however, the 
nonlinear correction is different. In fact, it is 270% times 
larger due to the Hall current having a dominant effect 
around the resonance. Moreover, the dispersion in the 
slow dissipative layer causes a further reduction in the 
coefficient of energy absorption, in comparison to the 
nonlinear regime alone. 

At the Alfven resonance dynamics can be described 
within the linear framework. Hence, using Eqs. (|59p and 
(l60l) we obtain that 



r Q = 



(r a + m) 2 + ' 



(63) 



Numerical verification of these results requires much 
more work than would first appear, and as such our next 
paper is to concentrates on this and further numerical 
analysis. 



VII. CONCLUSIONS 

In the present paper we have investigated (i) the 
effect of nonlinearity and dispersion on the interac- 
tion of fast magnetoacoustic (FMA) waves with a one- 
dimensional inhomogeneous magnetized plasma with 
strongly anisotropic transport processes in the slow dis- 
sipative layer (ii) the interaction of FMA waves with 
Alfven dissipative layers. The study is based on the non- 
linear theory of slow resonance in strongly anisotropic 
and dispersive plasmas developed by Clack and Bal- 
lai(author?) [21| and the theory of Alfven resonance 
developed by Clack et al. (author?) [221 ]. 

We have assumed that (i) the thickness of the slab con- 
taining the inhomogeneous plasma (Region II) is small 
in comparison with the wavelength of the incoming fast 
wave (i.e. kxo <C 1); and (ii) the nonlinearity in the dis- 
sipative layer is weak - the nonlinear term in the equa- 
tion describing the plasma motion in the slow dissipative 
layer can be considered as a perturbation and nonlinear- 
ity gives only a correction to the linear results. 

Applying a regular perturbation method, analytical so- 
lutions in the slow dissipative layer are obtained in the 
form of power expansions with respect to the nonlinearity 



parameter £. Our main results are the following: Nonlin- 
earity in the dissipative layer generates higher harmonic 
contributions to the outgoing (reflected) wave in addi- 
tion to the fundamental one. The dispersion does not 
alter this, however, the phase and amplitude of some of 
the higher harmonics are different from the standard non- 
linear counterpart (see discussions before) . Dispersion in 
the dissipative layer further decreases the coefficient of 
the wave energy absorption. The factor of alteration to 
the nonlinear correction of the coefficient of wave ab- 
sorption due to dispersion is 270%. Remember, however, 
that the nonlinear correction is multiplied by the small 
parameter £ 2 , so the effect to the overall coefficient of 
wave energy absorption is still small. 

Calculating the coefficient of wave absorption at the 
Alfven resonance confirms the linear theory of the past 
and verifies the approach taken to be correct. As our 
physical set-up of the problem (for the Alfven resonance) 
matches the typical conditions found in the solar corona, 
these results can be applied to it. The equilibrium state 
of the problem (for the slow resonance) can match con- 
ditions found in the upper chromosphere, where FMA 
waves may interact with slow dissipative layers, and if 
the reduction in the coefficient of wave energy absorp- 
tion persists to the strong nonlinear case (as with the 
long wavelength approximation found by Ruderman [27| ) 
dispersion may have further implications to the resonant 
absorption in the solar atmosphere. 

In a forthcoming paper, we shall theoretically and nu- 
merically investigate coupled resonances, which builds 
from the work in the present paper to obtain a more 
realistic model for a solar physical description. In the 
same paper we will numerically analyze the absorption 
of fast waves at the Alfven resonance as a possible sce- 
nario of the interaction of global fast waves (modelling 
EIT waves) and coronal loops. 
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APPENDIX: DETAILS FOR CALCULATION OF 
FOURTH ORDER APPROXIMATION 



It is straightforward, but longwinded, to derive qi 2 , so 
we skip all intermediate steps and give the result 



In the fourth order approximation Eq. (126p gives 



09 



-id 2 94 

" ae 2 



V 4 



dp. 



4 c 



dq 2 dq 2 
da 86 



d_ 

m 



p 0c v Ac \A\x dO 
1 2 \ dqi dq 3 dq 3 dqi 

2 q * (A1) 



We can rewrite the third term on the right-hand side of 
Eq. (|A1[) using our knowledge about the first three orders 
of approximation, so 



d_ 

00 



<7i<73 



2 

7^2 



(<?lQ 3 l +<7l933) 



2ik6 



+i(2q l q^ + ql)e^ k6 }. (A2) 



This equation contains terms proportional to e 2lke and 
e Atke , so we can anticipate the solution to Eq. (1A1I) to 



be of the form 54 = K (g 42 e 2l ' ct ' + (744c ) , where 34 
represents P4, W4 and 54. We calculate the fourth order 
approximation to demonstrate that nonlinearity and dis- 
persion in the dissipative layer generates overtones in the 
outgoing (reflected) fast wave. For brevity, we shall only 
derive the terms proportional to e 2%k6 , but for complete- 
ness we note that it can be shown that terms proportional 
to e 4lke are only present in the solution inside the dissi- 
pative layer. 

Using the continuity conditions at x — and x = xq 
we find the jump in the normal component of velocity 
across the dissipative layer to be 



W42 



iKiVA A2 



n e VA A2 



Pi (V 2 - v 2 Ai cos 2 a) p e (V 2 - v 2 Ae cos 2 a) 

- 2ikVA 2 @> / F~ x {x) dx 
Jo 

-2ikV(kx )& -^Hrdx. (A4) 
Jo F{x) 



-1 



<742 



(j — 2% 
V 16 (p, 



V A A A 



v*( Pe + A 1 )A 31 (i + n 2 ) 



Pocv\ c 



A\x 2p 0c v Ac \A\ 2 x 2 (<T 



Ax) \ Pe + A,\ 



96pg c < c |Ah4(<7-i)3( CT 




(A5) 



with f^42 = /(c), where f(cr) — > as a ~ > 00, is the 
contribution due to the Hall effect. As it is not essential 
for forthcoming calculations, its exact form is not given 
here. The substitution of <7 42 into Eq. (|2T|) yields 



2TrkV 5 A 4 



W42 



p 0c v\ c \A \ cos 2 a 

A.»£|A| 4 igcos 2 o 
Comparing Eqs. (|A4[) and (|A6|) we obtain that 



A42 = 1.279 x 



4 4 4 6 

V e T p cos D a 



n 3 V 3 k 3 x^ (p + ivf {p 2 + v 2 ) 



(A6) 



(A7) 



Here we have used the same estimations that were uti- 
lized for calculating A31 in the third order approximation 
and retain only the largest order terms with respect to 
kxQ. It is clear from this result that the outgoing wave 
becomes non-monochromatic in the fourth order approx- 
imation. We can also observe that the second harmonic 
appears in addition to the fundamental mode. 

This result parallels the results obtained by Ruder- 
man et a/.(author?) ^ andBallai et aZ.(author?) [3. 
However, Eq. (|A7|1 shows that the phase is inverted and 
the amplitude of the second harmonic is approximately 
30 times greater than theirs due to the presence of the 
Hall effect. Remember, though, that this amplitude is 
multiplied by a very small term, £ 3 , which means the 
overall correction is very small. 
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